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АВ5ТААСТ 


The digitizing procedure used at the Naval Postgraduate 
School was investigated for possible sources of noise and орке 
errors. Signals of known form were digitized through the 
Analog-to Digital Hybrid computer system (Ci 5000/XDS9300). 
Similar signals were generated by digital programs on the 
IBM 360/67. The resultant signals were analyzed by the com- 
puter programs UBCFTOR, which computed the Fourier coefficients 
of each block of data, and by UBCSCOR, which computed the 
Rewer Spectra of the signals. The power-spectral plots of 
the computer-generated signals were compared with the power- 
Bal plots of digitized signals. The analog-to-digitital 
process appeared to be relatively noise free. 

To further test the system, atmospheric а T апа 
wind velocity signals were digitized and analyzed under 
UBCFTOR and UBCSCOR. Plots of the time-varying spectra of 
these signals compared favorable with results obtained at 


vT digitizing facilities. 
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I. INTRODUCTION 


A. PROBLEM 

Random processes are often recorded in a fluctuating vol- 
tage in analog form. HO it is frequently more converne ne 
to analyze the signal digitally (on large digital computers). 
Bus, investigators are faced with the prob lene ei омеае 
а сопС1пцоц5 signal into discrete data samples, and with sub- 
sequent analysis of these digitized samples. There are 
several steps in the analog-to-digital conversion procedure, 
and in the digital analysis procedure, which allow. for errors 
Ea tor possible contamination of a signal with noise from 


external. sources. 


ЕЕ OBJECTIVE 

me Main goal of this study was to investigate digitization 
and analysis procedures for possible sources of noise which 
may be introduced to the true signals. iheroxeral isproer due 
used at the Naval Postgraduate School, from een of 
the actual signal to the computation of the Power-Spectral- 
Density (PSD), is quite complex and requires four computer 
Broerams. The next objective was to improve the routine 
procedures required in thís particular time series analysis 
technique. The final objective was to digitize actual | 
geophysical signals and compare the results with similar 
analyses of these signals undertaken at the University of 


British Columbia, by Boston in 1970 [Ref. 1]. ^ 





ТЕ ЕО 


A. DIGITAL REPRESENTATION OF CONTINUOUS, TIME-VARYING SIGNAL 

Random geophysical processes are often studied by TES 
cording a continually changing event as a continuous, fluctuat- 
e voltage, which corresponds linearly with the original 
ШЕСсесв. The actual geophysical variables cons ее ES 
in this study were small-scale fluctuations of air temperature, 
wind velocity, and time derivatives of both. 

1. Analog-to-Digital Conversion 

A randomly fluctuating voltage signal might look like 

the one in Figure l(a). In order to analyze the signal by 
ial techniques, discrete samples of the fluctuating voltage 
must be taken, and these are referred to as ведцеп аш ате л 
samples (vi). The requirement for sampling at equal intervals 
of time is set by the assumption, in most analyses, that they 
are equal time interval samples. The digitized samples would 
look like Figure 1(b). 


Б. Digital Sampling Theory 


implied Tin the digitization problem, E dle LION 


~ 


ИН well the sequence represents the original voltage. In 
Mining to Sampling Theory for the answer, three hypotheses 
eee made; | 

a. V(t) is a random variable defined for -»<st< о, 


b. The power spectral density of G(f) exists. 


15 ` 








a) Analog Signal 





/ b) Digital representation of analog signal 


Figure 1. Analog and Digital Signals 
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c. G(f) equals zero for frequencies equal to, or 
greater than B. B is the highest frequecy which can бе ге= 
seuved and f is defined as frequency. 

Further, if we let At be the sampling interval in seconds, 
so that Ats1/(2B) and if we have v(t) sampled at intervals of 
EB giving vi samples, i = INED f it can be shown 
that v(t) can be reconstructed uniquely. 36. if the power in 
v(t) is limited to a band less than V(2At)Hz, then, sampling 
at an interval At allows v(t) to be oe ппшде Буй 
Oof of this augument is given in [Ref. 2]. In this way, 
ШЕР requirement of sampling at least twice per cycle has not 
only been established but is entirely sufficient. | 
ES. Limitations on Frequency Resolution 
К е to real world limitations of not being able to 
collect a signal v(t) of infinite length, and because our 
actual signals are not necessarily band-limited, we have to 
medify many of our theoretical assumptions. Firstly, we 
те we can get a long record length, which is representative, 
at least over the range of frequencies we are СЕ ЕЯ а 
of the signal extending in Cime tO plus andi minus ini imiity, 
secondly, we can use "low-pass" filters to eliminate un- 
Memted high frequencies before the signal is sampled. 
a. Low Frequency Limitations 
If the signal V(t) contains low frequencies it 
will be very hard to distinquish then in the interval 0<Ғ<1/(2Т), 
according to Rayleigh's Criterion, [Ref. 2], where T is the В 


record length іп seconds. Іп this situation, when recording 


Џо 





a finite section of signal, we have, in fact, киел есы ше 
erpiginal signal. The recorded section oT даса шен ene sens 
the original signal. The term "block" refers to a further 
Beuncation of the signal into smaller sequential sections ven 
data, which, when added sequentially, will give a truncated, 
mp led, section of signal. Thus as T increases, or the 
length of signal examined ПО the lower will be the 
frequency which can be resolved. 
b. High Frequency Limitations 

The highest frequency we can resolve has been 
established as one-half the sampling rate. In other words, 
at least two samples per cycle are required. The high fre- 
quency limit is more commonly known as the "high frequency 
cut-off point," or just "cut-off frequency," (fo). Sometimes, 
Requirements of five samples per cycle are sev; however. chris 
added computational problem is in fact unnecessary and 
eeecrntially means that the high frequency limit of analysis 
Ши сгеазед by a factor of 2.5. If 1000 Ss were re- 
quired for a high frequency limit of 500 cycles, 2500 Samp/ 
Sec would ме in û high frequency lIImil OI I290 CVECleS. 
Usual sampling rates vary between 2 and 2.5 samples per cycle. 
This study, a sampling rate of two samples per cycle was 
Wa for the turbulence analysis. This rate had been 
established as satisfactory by Boston (1970), [Ref. 1]. 

C. Aliasing 

The requirement of simply sampling Леле 

frequency of interest at two samples per cycle is adequate 


ШІ the highest frequency of interest is, in fact, the 
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highest frequency present in the signal. When higher fre- 
mpemcies are present, the cut-off frequency becomes selec 
"Folding" or "Nyquist" frequency. The "folding" comes from 

the fact that higher frequency energy, "electrical energy" 

our case, is "folded" back into lower frequencies around 
mos point. To eliminate the problem it is necessary to either 
mle at higher rates, thus moving the folding frequency 
Emer, or to sharply filter the signal at the folding 


Brequency. 


В. FOURIER TRANSFORMATION = 

Many studies of random geophysical processes seek to 
@escribe the distribution of energy within a varying process 
memo unction of the frequency of fluctuation of the variable 
being GULES: sometimes when analysing turbulence, it is 
Med to determine the distributation of energy in the 
Mxdly fluctuating velocity as a function of frequency. The 
problem is, thus, basically one of transforming data from 


the time-domain into the frequency-domain. 


ШИ Еспгіет Transformation of a Continuous Signal: 
w ourciler Integral 


One of the most often used ies for this trans- 
formation, from the time-domain to the frequency-domain, is 
eh the use of the Fourier Intergral Transform. Accord- 
ing to this procedure, a periodic sinusoidal signal (shown in 
analog form in Figure 2a) is transformed from the time-domain 
into the frequency-domain through the transform-function 


-janft 


V(f)= fv(t)e dt (1) 


Where v(t) = E Sin 21 f¿t and j=v-1 


ЧК - 





as AMPLIFUTE 











Fourier Coefficient 
Fourier Amplitude 


Transform 





‚Time Г Lo Frequency 
“E E 
LP 
а) Periodic Sinosoidal Signal b) Fourier coefficient 
| Of kransiormee 
signal 
AMPLITUDE 
VoLTS 
f, fa fa fa 
frequency 
c) Complex periodic Signal d) Fourier coefficients 
| ОҒ transformed 


signal 


Figure 2. Fourier Transformations 
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This transformation is shown graphically in Figure 2b). A 
Benplex, almost periodic signal eomposedzofizeseyer az are 
waves of differing amplitudes and differing frequencies would 
be similarly transformed into the frequency-domain, as shown 
in Figure 2d). | 
2. Fourier Transform of Discrete Data Signal 
When dealing with digitized signals or discrete data 


the finite form of the Fourier-Transform must be used: 


N-1 -12тҒ14А% 
Ук) = ВЕ ме —e«f«o (72) 
1=0 í 
where v is a complex variable and K=0,...,N/2. For real 
data vi where i=0,...,N-1, the sine and cosine transforms 
becomes: 
- N-1 
a(fy) = Att Va Sin 21f;At | (За) 
1=0 
N-1 T» 
ВЯ = ш уу Cos 2mfj4t (3b) 


erence 2 suggests that, computationally, this would follow 
a flow diagram as in Figure 3. The value of f is computed as 
 Шпсесіоп ог the number of data points and the Tene Conon 

the récord in seconds (T). In other words, increasing the 
шлога length decreases the frequency difference between co- 
cients; or a longer record will give eT estimates 
closer together in frequency. Making the substitution for 


f the sine transformation becomes: 


K? 
N-1 A iK 
a(KAf) = F У; Sin 27— (4) 
1=0 N 


This approach to computing Fourier coefficients has been 


limited in the past, due to the computational time required. 
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Compute 
^X Cos O 
Vai Sing 









Accumulate 
Both 
Terms 








Figure 3. Flow Chart Showing Fourier Transform Procedure 


for use with Discrete Data 
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The number of calculations required Increases OS ие 
of the number of data points increases. The integral trans- 
form, in the past, was limited to studying theoretical 
functions, such as sine waves and square waves, which were 
relatively well-behaved functions. For cases such as these, 
even the discrete transform could be ned to’ readily Obi; aan 
res Pourier-transform of a signal. With the introduction of 
the Fast Fourier Transform (FFT), algorithm, computation time 
E been greatly reduced. 

ШЕШ Fourier Transform of Truncated Continuous wave Form 

If a signal v(t) exists only for the time interval 

from 0 to T seconds, and is zero at all other times, its 


Fourier-transform is: 


T -jamft 
V(f) = f v(t)e dt 2254 
С | 


If v(t) is repeated in intervals of period T seconds, the 
frequency difference between coefficients will be 1/T Hz. The 


Kth coefficient will be: 


T -j2nKt 
Ук =i / v(t) e 1 dt (ба) 
TIG 
(ог the discrete case: 
К: 
N/2 -jen—i , 
ШЕ - 277 у е Е (6b ) 
К Т о з=] | 


ШЕ Convolution of Continuous Signals 


Ihe finite-transterm of a finite-length tcine series 
can be viewed as the product of a finite-length rectangular 


function C, _(t), times an infinitely long-time history v(t). 


72 
The finite transform of ví(t) becomes: 
T -janft 
E ые dt (7) 
0 ге | 
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Си псе products transorm into convolution Teee conyoluted 


Fourier-transform becomes 
Vo(f) = VW(f)*Gpyo(f) (8) 


where 
T -j2nft 

@т/2 = гды e dt (9) 
Figure 4(c) shows the Fourier-transformation of the rectangular 
function and Figure 4(d) shows the convolution of a sine wave 
Ж ЕГІ phe rectangular function. This convolution problem 
arises when switches are opened and closed while recording 
the data. The side lobes of the convolved function can be 


minimized by allowing the time length of the record to be 


sufficiently large. This decreases the interval O to 1/T. 


C. POWER SPECTRAL DENSITY EUNCTION 
in Methods of Computing Power Spectral Density 
There are three methods which may be used to compute 
mhne power spectral density of a signal. They are: 
a. Direct Fourier Transform Method 
The Fourier transform of the signal is computed 
and from this the mean value squared is determined. 
b. Analog or Bandpass Filter Method 
The signal is put through a bank of bandpass 
filters and each filter output is squared and integrated. 
Cr Auto-correlation Huneticnelletchod 
The Auto-correlation function of the time Series 


is computed, and then its Fourier-transform is computed. 
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a) Rectangular function 





b) Truncated sine function 





-£o 


d) Fourier transform of truncated sine function , 


Figure 4. Effect of Truncation of Sine Function 
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Historically, the direct method Was develop eae ee 
but could be used only on nicely bebayed theoretical signals. 
When less well behaved signals where analyzed, smooth 
(deterministic) theoretical functions were replaced by discrete 
data points representing the signal. This necessitated using 
meee discrete Fourier-transform; however, for large data-sets, 
е vime for calculating the Fourier-transform was prohibitive. 
Ше computer program FTOR utilized a variation of the direct 
method employing the Fast Fourier Transform (FFT). | 

Using the sampling rates 1000-5000 samp/sec, (SPS), 

a vine the ability to select the block sample length and 
meemmumber of blocks desired for a particular run, no pro- 
blems were encountered in which filler data, usually Zeros, 
паа to be inserted into a. block. The block size was always 
Seen as ah integral power of two. 

2. Typical Power Spectral Density Functions 

Ере Tower Spectral Density(PsD) plot for ralidom dala 
shows the distribution of electrical power within the signal 
as a function of frequency. Several characteristic PSD plots 
are encountered. Figure 5(a) is a PSD plot of.a pure sine 
wave. It gives the Dirac-delta function which implies that 
che power at the sine TA is infinitely large, and жеге 
IL other frequencies.. Figure 5(a) is a PSD of a Gaussian | 
Eom signal. The PSD of this Signal is constant. Figure 
5(c) shows the PSD of a random signal carrying a sine function 
Smet. The PSD plot is the sum of the power spectra of the. 


random signal and sine figured separately. 
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PSD (vots na) 


to Frequency (8 2) 


a)_ Sine Wavê 


85 (мо УМ 2) 


Eyrequeney (Ha) 


b) Gaussian Noise 


POD (Volts Az) 


Lo Frequency (Kz) 
c) Gaussian Signal Plus Sine Wave 


Бириге 5. Characteristic PSD Plots 





3. Problem of Sina lG Sided E 


The Fourier coefficients under the transformation 
-janft' 


OD Za y(t)e dt (10) 
Beet for both positive and negative frequencies. the co- 
efficients can be transferred to a single-sided frequency plot 
С Emp ly doubling each coefficient aS iL 15 plotted only a- 
long the positive frequency axis. If S(f) represents the two- 
sided spectral density function and G(f) represents the single- 
sided spectral density function, the single-sided spectral 
density function equals twice the double-sided density function: 
G(f)=2S(f). Likewise, if a watt meter was used to find the 
Beer in a signal as a function of frequency, the values ob- 
tained would have to be divided by two before plotting оп а 
а-а spectrum. 
4. Power and Energy Signals 
The average : electrical power in a fluctuating 


voltage signal v(t) is given by 


P = lim 1f? |(t)]? at | (11) 
аа 


Ref. 3] defined a power signal as опе which is FoF all 
practical purposes, infinitely continuous. Energy signals 


are pulse-like in form and are given by: 


г 


Е = Л) |v(t)]^ at | (12) 
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The power in a periodic signal is shown from 


Parseval's Theorem to be: 
t'+T 


p lv (t)| f at 


1% 
Thus, the power in a continuous sine function, v(t)=sin 
would be ; 
2 = 
poc Lo Sasso upE disce ле га 
ТІ 


If a real wide band signal were to be analyzed by using 
able filter of band width of óf, a plot of power versus 


frequency as shown in Figure 6 might result. Here, the 


(13) 


Se. 


(14) 


a tun- 


real 


weet from zero to an upper frequency was determined, divided 


by two and the values folded back into negative frequencies. 


5. Power from the Fourier Coefficients 


Parseval's Theorem leads to the definition of power 


in terms of the Fourier coefficients 


DE а> 


да 
пи --= УК 


XP 8 


Thus, knowing v(t), the average power can be computed . 


average power in a one volt, 60Hz sine wave is: 
GE E O COE 


Р = 


T 
Lo 


зе 


Sin? 377t dt = 1/2 watts 


Bez the Fourier coefficients: 


| 
O 


У_1 = 1/2 Уу = 1/2 Ук„у = 


1/4 + 1/4 
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Figure 6. Two-Sided Power Spectrum of Wide Band Signal 
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Figure 7. Power Spectrum of 1 volt 60 Hz Sine Signal 
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If a watt meter was used to check the power of this 
sine signal, one-half watt would be found at 60Hz, and zero 
Watts at other frequencies. If a plot of power yersus fre- 
quency were made, one similar to Figure 7 would resulte s 
ШЕ values have been plotted for positive and negative fre- 
quencies, the total power has been iced by one-half and 


the values have been plotted. 


D. FAST FOURIER TRANSFORM | 
|. Computational Economy Afforded. by FFT 

The Fast Fourier Transform (FFT) algorithm is a 
ШЕ песпос for computing the finite Fourier transform of a 
Bes of N data points. The FFT computation requires N 
log, N computations, which leads to a substantial saving in 
Computer time over the old method which required N* Operations 
the FFT economy of computational effort is greatest for large 
pues of N. 

Reference [4] gives the general computational 
femme for figuring the FFT. It is pointed out that the 
PFT is general in that N is not necessarily a power of two; 
ver, Dy Selecting N tO be a power or two, further computa- 
onal savings result. When this requirement is met, the 
FFT algorithm is essentially a successive doubling operation. 
is for N=2d , only Nj multiply-add operations would be Dez | 


quired under the FFT assuming the necessary complex exponential 


table of values has been computed in advance. 
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2. Importance of FET in Turiu Analysis 
Computationally, the FFT 15 ап іпіісвренсарле атана 


the PSD investigation of turbulence: Due to the extremely 
Шел sampling rates necessary Tor studying high wave number 
processes, enormous amounts of data are collected. A five 
minute section of a fluctuating temperature signal sampled 
at a rate of 4000 samples/second (SPS) generated almost a 
million and a half data samples. Using the old method for 
muting the Fourier transform for this quantity of data, 
E thousand billion operations would be required. Clearly 
Me” time element in this procedure would prohibit the 
 ШопТасіоп on any but the fastest computers. Turning to the 
FFT the operation would only take about 21 million operations 
or an improvement in excess of 5 orders of magnitude in 
eemputational time. The PPE Gus (brings CMe SEMT of 
Mmemepulence Signals within the bounds of computational 
Feasibility. 

The ability to determine the power spectral density 
function which is computed from the Fourier-transform and 
which is such an important aspect or Да ле ие 


made possible through the use of the FFT algorithm. 
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ТТТ. NAVAL POSTGRADUATE SCHOOL DIGITIZATION FACILITY AND 


SPECTRAL ANALYSIS PROGRAMS 


A. NPS COMPUTER FACILITIES FOR DIGITIZATION AND PSD ANALYSIS 

The Naval Postgraduate School has two sophisticated computer 
systems which are used in the digitization and PSD analysis 
procedures. One is a Hybrid system which consists of an 
analog computer, COMCOR Ci 5000, which is electrically inter- 
faced with a digital computer, XDS 9300. The XDS 9300 con- 
tains 34K bytes of core storage and uses an octal number base. 
Hes base consists of binary numbers made up of 3-bit digits. 
The second system is an IBM 360/67 which has a core storage of 
762K bytes. It uses a hexidecimal number base which consists 
of binary numbers of 4-bit digits. 

i Hybrid Computer 

аск пи compter facility is located om thew: th 

floor of Spangel Hall. The equipment disposition is shown 
in Figure 8. Though the individual user performs all computer 
operations, technical assistance is availabe from ene 
Electrical Engineering Computer Laboratory staff. This 
facility is used for the analog-to-digital- conversion of 
senals. 

2. IBM 360/67 Computer 

The large, digital computation facility S oad on 


wei ground floor of Ingersol Hall. Computer Center personnel 


handle all computer operations. Individual programs are input 


æ 


E 





Катттови лодзпашо5 ОАО Н 





pe 





to the computer through a card reader andi regust E o 
subsequent input of programs, tapes, disk memory, etc. are 
submitted "over the counter" to computer personnel. A 24- 
hour operating schedule is usually maintained. All PSD 


Ensvses of digital tapes were performed on this computer: 


SINE YBRID COMPUTER SYSTEM | _ 

Computer time on the Hybrid computer system is signed for 
in advance in the Electrical Engineering Computer Laboratory. 
The whole system is operated on a self-service basis; however, 
mmemecomputer center staff is available to assist during 
working hours. : Computer time was easiest to obtain early 
in the term or between terms when the work-load was entes 
The facility is available 24 hours a day and as one becomes 
more proficient in equipment operations, night or week-end 
operations can be scheduled if week-day operations are precluded. 

1. The Analog Computer (Ci 5000) 

The Analog Computer contains the input points for 
raw Stenals input. It also has control features for signal 
tification, selection of sampling rate and starting and 
stopping of the digitizing procedure. The two removable 
топ boards are the heart of the system. Individual 
patchboards have been reserved solely for analog-to-digital 


conversion (board #24). 


33 





раводАеу pu? сраводцојед 0006 





ЈЕ 


J 


е 


6 


5 


Эли E 


Озуму 


34 






цозтмс чо 1IMOJ PUB PABOQAYSA 0006 


A 


b 


О 








The logic board, the smaller of the two равеспран ин 
was utilized to. set the sampling rate of Une anevocece ae 
ihe logic board with its patch is Shown an Hicures a = 
The sampling rate was changed by turning a counter-like number 
=m play, below and to the left of the logic board, as shown 
ЕБР риге 17.8 Тһе counter operated a voltage divider which 
mes clectrically connected to a resistance input on the logic 
board. The resistance value at this point was changed by 
Meoerling a resister in one of 4 positions. The resistance 
value was then divided By the counter setting plus one o 
give the sampling rate (in samples per second) TT 

Example: 100ke _ 4Kc or 4000 samples per second. 

(24+1) 
No external leads are connected to the logic board. 

. The analog patchboard is the input point for all 
ekternal signals entering the analog computer. The input 
Signal from a tape recorder or signal generator, aA 
through signal conditioning equipment is input to one of the 
analog board amplifiers. (Figure 15 shows inputs going into 
amplifier A001). Various gain factors can be applied to the 
ut signal to bring it up to a value optimum for НО 
of ne analog computers dynamic range. Figure 10 shows the 
сага ог the analog computer which is used to control the 
enred operational mode of the analog computer. 

2; Digital Computer (XSD 9300) 

The heart of tne digital computers Sens CIC тох 

Data System 9300 Central Те. Unit (CPU buc ле туе 


units, line printer, card reader and teletype unit are inter- 


faced with the CPU. Figure 11 shows, diagrammatically, the. 
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Bquipment used in digitization. The каре ету сы ern 
Ша Гог data input from tape to computer or гог аара е 
SM from computer to tape; the line printer can be used ier 
ШО ат listings and printing results; the card reader as pro- 
ЕСІП input for compiling programs; the teletype unit ror input 
of parameters necessary to control the "multi-channel analog- 
to-digital Program." The CONTROL CONSOLE of the XDS 9300 is 
Mco compile and to initialize the program. After the pro: 
em has been EN спе teletype 1s cu select one of 
ВЕ: Gigitizing options. The digital and analog computers are 
connected through the Analog-to-Digital Converter (ADC), (Fig.12). 
3. Multi-Channel Analog-to-Digital Program 

several programs have been developed by the computer 
шерогабогу staff to control the XDS 9300 system during a 
nmel digitizing operations. A card deck of the latest re- 
Vision of the program is maintained in the computer lab. This 
meram can be input from either cards or from a tape which 
has nad the machine language program stored on it. The com- 
Ме Sequence takes about five minutes using the card input 
and only about 30 seconds using the tape input. 

enee Che р опат нас БЕС come о EMIT 
Program Control Options can be selected by typing one of the 
options followed by carriage RETURN (C/R) on the teletype 
(see Figure 13.). The Program. Control Options are: 

1. Enter new parameters. (NSAMP, NCHAN, NREC, ITAPE) 

2. Start dlgitizing the апалов отары mare . 
wit ization actually starts when manual switch БӘТ оп od PN 


ПЕ оона to up position). 


E) 





ANALOG DIGITAL 
TAPE TAPE 


n ( HILTER С:5000 М ey 9300 
UNJT 





Figure 12. Block Diagram of Analog-to-Digital Conversion 
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3. Write End-of-File on tape. 

Ц. Rewind tape to load point. 

5. Skip files (the number of files to be skipped 
јео Гједа by teletype inputs of a four digit nurni e ЕС SK 
five files 0005 as shown in Figure 14). 

6. Print digitized data (the number of lines to be 
printed specified by teletype input of a four digit number, 
BEEN rn. ten lines of digitized дага 0010, 

T. Actuate the Digital-to-Analog subroutine for the 
Teme single block of data encountered and input CHORDS into 
the strip chart recorder (a subsequent computer comment types 
"start strip recorder and type C/R - carriage return"). The 
strip recorder must first be connected to the digital-to-analog 
output on the analog patchboard on the Ci 5000. 

After the multi-channel digitization program has been 
compiled and the input light on the teletype is on, tape 
parameters input through the teletype are: 

МВЕС - © Sets the maximum number of records to be 

digitized (see Figure 14). | 

^ NSAMP - sets the number of digitized samples contained 
in each record on the seven-track tape 
(see Figure 14). 

ITAPE - sets the magnetic tape unit into which the 
digitized samples are sent. This number must 
be the same as the dial number selected on 
the front of the tape drive unit (see Figure 14). 

NCHAN - sets the number of analog channels сс бе 


digitized (see Figure 14). 
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С. ANALOG- TO-DIGITAL OPERATIONS 
1. A ESA 
а. Analog Tape Recorder 

The most commonly used record/playback device in 
the Oceanography Department is the Sangamo tape recorder 
Model 3562, which can record up to 14 tracks of data, and 
two edge-tracks of voice on one inch magnetic tape. Two modes 
EE Deration are possible: direct and FM. It is general 
КОШ се to only use the FM mode, which records from D.C. to 
some upper limit depending on the tape speed. The FM mode 15 
ЕЙ ПСІЗе Ігев. Direct electronics are used for higher 
meeauency data when no D.C. information is required. 

The magnectic heads and tape rollers of the tape 
recorder should be стеанеч Wilh awseopropyl alconolmafter pro- 
longed use to reduce noise during playback. 

De ЕНИ етв 

KHRON-HITE, Model 3321 filters have been most 
widely used in the past for filtering out unwanted high and 
Шгепцепсіесв from a signal. They can either by connected 
Singlely in a low-pass mode or in series as a band pass fil- 
ter the filter mode has been selected the raw signal 
ШАШаприс to the filter and the filter output connected to the 
Ос patchboard as an input to the analog computer. 

c. Analog Patchboard 

The analog board should be patched as shown in 
meure 15 if two channels are to be digitized. If only one 
channel is desired the signal inpite. should be aioe. 


amplifier AOOQl only. Gain factors may be patched as desired. 
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The outputs of each amplifier во горе е 
oscilloscope input points and the yc points 500 and 
7501). If an analog playback of the digitized signal is 
desired, inputs from T420 and T421 (see Figure 15) must be 
made into the recorder inputs (1-8). 

di. Logic Patenboard | Ў 

The patch for the logic board is shown in Figure 
Nils 15 for continuous mode digitizing in which digitization 
Mota rted by throwing switch Dsl. The only logie setting 
eee ary is the selection of the required resistance value 
(Figure 17) which in conjuction with the wheel counter sets 
Sampling rate. 

О Обе Тозе је 

The "scale illumination" dial energizes the oscil- 
loscope (see Figure 16). The sampling pulse was ШЕШЕН 
patehed into "channel 1" by throwing switch Y1 to the left. 
This permitted continuous monitering of the sample pulse fre- 
quency and width. Dial DFOO (see Figure 17) should be set to 
Ms. 1 and the width of the sample pulse adjusted шоа је аза. 

EE Energizing the Ci 5000 Computer 
If the power on the CI 5000 has been secured (see 
power light in Figure 9b) it can be turned on by depressing 
the "on" switch. Next, the buttons KEYBOARD, POTSET and 
RESET are depressed in that order (see Figure 10). 
3. Energizing the XDS IO 
The step by step procedure below should be followed 
in energizing the Hybrid computer system. The latest update 


to this procedure is kept on the XDS 9300 Control Console. 
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а. XDS 9300 Power on (if power is off). see Figure 18 

1) Depress IDLE 

2) With RESET depressed press POWER 

3) Turn teletype switch to ON 

4) Press CLEAR and CLEAR FLAGS simultaneously 

b. Energizing the line Printer (if the READY light isn't on) 

1) Insure that POWER button on JM printer is on. 
ШИ ПОС, press POWER button. 

2) When lower half of POWER is lit (within one min. 
after POWER pressed) press READY, NOTE: READY must be off to ad- 
vance paper. A mote page is advanced with TOP of FORM and the 
paper is advanced a single line with SINGLE SPACE. 

с. Card Reader (if A/D program being input from cards) 

1) Place BOOT card in front of the JOB card in the 
A/D deck. 

2) Put cards into card reader (see Figure 19) С 
down, top outward (toward you) so that column 1 is first into reader. 

3) Put card reader press in order: POWER and 
START. 

d.. Program Compile and Load 

1) Ready the Line printer (sequence b above). 

2) With cards in place ready the card reader 
(sequence c above). | 

3) At 9300' console turn off (by depressing) all 
SENSE switches which are lit. | 

4) Press in order: ЕРТЕ: RESCH ЕНЕ ARDS. 

5) A JOB is printed out by the teletype indicat- 


ing compiling has begun. 
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6) Upon compilation, the teletype will type 
END OF ASSEMBLY. Shortly, thereafter, the INPUT light will 
n, indicating the tape parameters can be typed im 

e. Equipment NOT-READY Condition 

1) If a device is not ready when needed, it will 
be called out by a teletype message. Simply making the device 
ready will clear the pause in the operation. 

2) In the event of a FEED-CHECK failure, inspect 
the card that caused the problem (the one on the bottom of the 
mut hopper).. If the card is found to be damaged on the 
leading edge, repunch the card and replace it. Put the un- 
read cards back in the hopper. Ready the card reader by press- 
ШЕ ГБЕСЕТ апа START. 

3) Line printer not-ready may becaused by its 
Demise out of paper. | 


f. XDS 9300 POWER OFF (after normal working-hours 


only ) 
1) Press IDLE. 
2) _ Turn Teletype SWILElN to OFF. 
3) Press IDLE. | 


4) With RESET depressed, press POWER. 
Ц. Mounting Magnetic Tapes | 
It is very important to use currently certified tapes 
Weer are free of БЕТ элеп errors. oCreawemes, tears, Сш 
spots, etc. cause tape-writing errors. The only recourse is to 
write an END OF FILE on the tape at the end of the bad file 


пик commence digitizing on the next Section Oi eave. 
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7-Track Magnetic Tape Drive Unit 
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Figure 20. 





Figure 20 shows one of the two tape drive units which 
па Бе selected for tape mounting еони оте о ет ние 
MANUAL READ. Before mounting the tape a circular plastic file 
WRITE RING must be inserted into the circular slot on the back 
of the seven-track tape being mounted. Mount the tape. If 
the ring is in place the FILE PROTECT, Figure 21, light on the 
ШЕР том of indicator lights will go out, indicating that data 
meee written on the tape. Tape threading instructions are 
Ма on the inside of the protective doors enclosing the 
tape reels. Once the tape is secure around the take-up reel, 
should be wound around at least six more times ЕР 9 
vm lape from Slipping off the take-up reel, when the main se- 
ле latch is closed. With the securing latch in place and 
the tape-idler-arm down, depress FORWARD-DRIVE. The tape 
Will move forward and stop unen iC redches tiem ime bal eE р. 
ШИШ ОО the LOAD POINT. The indicator light LOAD POINT on the 
ES Tov of indicator lights will go on. After the tape has 
been mounted and advanced to the LOAD POINT, the tape density 
is selected by turning the DENSITY dial to either 256 or 556. 
This determines che recording density iN Dits pel теа 
normally set to 556. The TAPE UNIT dial is turned колу опе 
ОГ several numbers, 1.е. 1. (The teletype input parameter 
ITAPE must match this number, i.e. ITAPE = 1.) Finally, the 
mode dial is turned to automatic. The TAPE DRIVE UNIT is 


раду. 
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Figure 21 {-Track Tape operating Dials and Buttons 





5. Variable Tape Digitizing Parameters 


Опсе the energizing and sampling sequences пате сеје 
completed (the INPUT light on the teletype should be on), the 
merece options are entered through the teletype Torey ne 
Sangle digit number and then depressing the CARRIAGE RETURN, 
C/R button. In order to input new parameters, "1 C/R" 
moma be typed. | 


Then the following example parameters could be typed 


in: 
NSAMP = 2048 (C/R) | 
МВЕС = 100 (C/R) 
МСНАМ = 2, ITAPE = 1, NDEL - 2000 (с /Н) 


mé computer would then type OPTION (Il), soliciting a further 
response from the user. If analog-to-digital collection was 
GO commence, the programmer would type "2 (C/R)." If the 
oes computer switch DSO was in the UP position, digitization 
would begin when DSI was thrown to the UP position. 

As would be expected, a high sampling rate would fill 
up 100 records faster than a lower sampling rate. Thus, if 
a sampling па ено и SES Mad Deen selected ава ие аселе 
баре parameters had been selected, the digitization would 
result in 204,800 samples being collected and a total signal 
length of 102.4 seconds. 


2048 SAMPLES 


RECORD * 100 RECORDS = 204,800. SAMPLES 


204,800 SAMPLES ` | 
5000 SAMPLES/SEC. x 2 CHAN 7 21.2 SECONDS 
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There would only be 102,400 samples of each channel; however, 
maere were two channels being digitized simultaneously) > US 
E Knows how digitized data for one and two channel diede 


would be formatted on a seven-track tape. 


D. CONVERT PROGRAM 

The CONVERT PROGRAM tonverts the seyen-track octal base 
data samples into nine-track hexidecimal samples. The basic 
Шиигассег ог the data samples being in blocks on the magnetic 
раре is maintained; however, two additional parameters (see 
Figure 23) are affixed to the beginning of each block by the 
CONVERT program. These are the numerical values KMAX which 
is set equal to the maximum number of samples per block and 
NCHAN, the number of channels of analog data digitized on the 
seven-track tape. This change in the block format, the change 
in number base and the addition of two new parameters is 
Meeeosary to put the tape in the proper format for input into 
the FTOR program. 

Figure 24 shows the CONVERT procedure in flow-chart form. 
Ses program processes all of the records ina file and the 
number of files to be converted is specified in the IF state- 
ИЕ which compares the number of files completed were 
number specified for processing. The number in parenthesis 
is set to the number of files to be converted. The JOB CONTROL 
LANGUAGE (JCL) cards which follow the CONVERT PROGRAM specify 
Aten files on a multi-file CHOC are to be processed. Jones 
[Ref. 5 Appendix D] gives the typical JCL cards Used voor. 


vert five files. 


m 


~ 





All data from 

single channel 
digitization is 
storred secuentially 


single one sample following 
channel WM the next, | | 
input " я 
QA iim i | oe 
| E. 3 RECORDS кој 
| Flle ONE FILE File 


~ 


a) Single Channel Tape Formatting - 3 Records with 8 samples 
per Record (N=23) 


Data from two channels 
digitized simultaneousl 


from alternate channel A N 


Number represents channe 
two number of data sample. 


Channel 
ШЕ TAA + 


input Sm! 
| 3 RECORDS 


ONE FILE 

















, 


b) Dual Channel Tape Formatting - Records with 8 samples per 
Шесога 


Figure 22. 7-Track Tape Formatting for Single and Dual 
| Channel Digitization 
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9-TRACK HEXIDECIMAL ТАРЕ 


Figure 23. 9-Track Tape Formatting for FTOR Program 
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DIMENSIC OS. IDA TZ јава тон 
ВЕ м1 м0 2 | 
REWIND 4 

IEND=0 
FACTOR=100.0/(2**23) 
(НЕСІ-2746 

NCHANz 4 

X«MAX=2J428 

J=0 









DS END= 40 EROS) DAT 
3| rormat(d2(6144)) 


ы / 7 


CALL FORM(ADAT PaEcL) 


SA TERS 


J 
10X,RECORD М0.- ,15) 


E 


10X,RECORD NO.= ,I5) 





DIR ARIAS | | 
4] DATCYJATOAT(I)* FACTOR Б 
7) 


Ce To URITE (6,66) (dat (1), 124 , 2048) 
N 66 РОНМАТ(1Х,8Е16.8) 


| МЯДТЕ(4) КМАХ, МСНАМ, АТ 

I NER 
CIEN 

E 
| »51 FORMAT(*0“,SX,READ ERROR,RECORD NO.=',1S) 
CO TO 31 | | 

ШАІТЕ (6,41) Ј--»>41 РОНХАТ(“0",5Х,ЕМО OF TAPE,RECORD NO.=',15) 
END EILE 4 

IENC- IEND«1 


IEND<NO, 
ЕЙСЕС 8 
PROCESS 












OF 







IEND«NO.OF í 
FILES TO 
PROCESS i 


REWIND 4 





READ £4) KIAX ,NCHAN,DAT 





м817:(5,82)/0АТ(1),1-1,2248)-1-82 FOR™AT(1X,8616.8) 
STOP | 
END 

ICU Cards Follow 


el’ For files 72272 тас Tape 
ӘСІ For Files To Зе Written 0n 4,-.%-.Tracr Tape 
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Figure 24. CONVERT Flow Chart 





Е. NAVAL POSTGRADUATE SCHOOL PSD COMPUTER PROGRAMS 

Seyeral computer subroutines are ‘available on the IBM 
360/67 for computing the Fourier transform of digitized in- 
put signals. Most programs, however, are designed for the 
put of relatively small data sets, usually from Hollerith 
cards. These programs could be designed to compute the 
Fourier transform of large amounts of data on magnetic tapes; 
however, a very powerful series of programs is available in 
the Naval Postgraduate School computer library which will 
compute the FFT.of large data sets. These programs are UBCFTOR 
UBCSCOR and UBCFCPL which are stored under files one, two and 
three respectively on tape NPS 216. The only particular format 
requirement for input data is that the digitized samples must 
be arranged in groups of' an integral power of two samples in 
each "record." A ee 2d samples where j is an integer.. 
Many records may be found on a single magnetic tape (see 
Figure 22). Since this FFT program package is especially 
well suited for analyzing large digitized data sets stored on 
magnetic tape, it serves as an essential component of the 
Naval Postgraduate School's signal processo capability. 
Figure 25 shows schematically the complete ion and 
PSD analysis procedure using these programs. 

The following description of the FFT package was taken 
from Dobson [Ref. 6]. Editorial changes were made to update 
the information with the Naval Postgraduate School facility. 
The programs were written by 2. Е. Garett and JI. Re Wilson 


while students at the Institute of Oceanography at the 
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migure 25. Schematic Diagram of Complete Рава лане 
| and PSD Analysis Procedure 
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Шаіуегсісу о” British Columbia. The UBC-IOR опанке 1m 
written in Fortran IV, instructs the IBM 360 to read the 
ital data tapes, store a block of data into the computer 
memory, compute the FFT of the data, and then store the 
resulting Fourier coefficients on another tape. From these 
coefficients the UBCSCOR program computes the spectral density 
от single SN or the eross-spectral density between 
selected pairs of MS Ihe spectral density values may 
be Pont put" in either tabular form or in graphic form on the 
Calcomp plotter. The UBCFCPLT plots the amplitudes of the 
rier coefficients as a function of frequency. 
ШО Pourier Coefficient Program UBCFTOR 

This program calls several subroutines which read 
in data from the digital tapes generated in the A/D phase. 
The maximum number of ar Points ато сап рел соге ояе 
computer memory under this program аге 8192 or 213 samples. 
NS restriction was imposed by PKFORT, the subroutine which 
EE Uutes the FFT. The details of this subroutine are listed 
under PKFORT in the Computer library. TOI OpLIMAL EIGEN CNC 
merORT has been designed to read in data blocks which con- 
tain an integral power of two (i.e. 24) senec PKEORTETS 
then called and the Fourier transform is computed, which re- 
IES in 24-1 Fourier coefficients. Since опту ollie сю 
data is analyzed at a time, the original signal is divided 
into ort, sequential sections of signal and a new length 
E гесога Т В The resulting coefficients are stored 


on another magnetic tape and they become the basis for the 


OS 





Epeetra. For this block of data, the hicho rE neo 

meme which a coefficient has been computed Ts еиен quite. 
Eueney fs/2, where fs is the sampling tuos pd The lowest 
Frequency which can be resolved for this data block is 1/T 

or 2d/fs. 

The Fourier coefficients are written трон се ке спин 
“coefficient" tape which also contains identifying information 
Such as the block inner, sampling frequency. nunberzarzscmeik 
Beach block, etc. Upon completion of the transformation and 
writing sequence TOF Ole block, the next block of data ls ресе 
Bon спе digitized tape. The sequence is repeated for as many 
mels aS specified on the input data cards, until an end-of- 
mere mark is sensed on the digitized tape or until a blank 
card is.encountered in the input data cards. 

2. Spectral Analysis Program - SCOR' 

Onee the Fourttemecoctiiicienvus паме been computed, 

Mie next step is to convert these values into spectral values. 
The program SCOR reads the coefficient tape generated by FTOR 
E averages the PSD values over 32 bandwidths. The data 

сага section of the SCOR program ЕрЕсл РЛЕБИ CRC UME OL 
channels, the number of blocks to be analyzed within a file, 
еп оЈоск number where the analysis is to begin, whether the 
Spectra for a single channel or the cross-spectra between 
channels is to be computed, the type of bandwidth desired 
(constant or logarithmic), etc. Just as with FTOB procedures, 
SCOR reads and analyzes E block of coefficients at a 


time. 
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If smoothing of the coefficients is deziro PE ee 


Enmoothing functions may be selecueda 


The "hanndug' opti 


performs a three point running average on the data with weights 


of -1/4, 1/2, -1/4. References 2 and 6, contain further de- 


pus of the smoothing functions. 


If individual Fourier coefficients. of two channels are 


By + il, and E, + E where i= 


2 р 


(-1)1/* 


and t equals the block 


length in seconds ( and 6F=1/t where Sf is the bandwidth between 


Fourier coefficients) then the power spectral density is given 


By: 
т(в + iia) 
a (E = e 
: 2 
2 2 
т(В- + 15) 
апа dro (1) ( 2 2 
2 


2 2 
Ry t a 


eof 


2 2 
Во + 15 


eof 


The cross-spectral values are given by the co-spectrum: 


q A (ВВ + ТТ.) 1208562 
0 2 26Р 
ema the quad-spectrum: 
БОЕ = ped Sow o -m* 
T(RR, 211) a 271 
ut 21 21 
2 28f 


, 


Mier tactor of 2 in the above equations makes the integral 


Mater the power spectrum over positive frequencies equal to 


Em Signal variance. 


Tase corrections can be made at this 


Meee Со correct for instrument phase shifting and then re— 


Euculating the cross-spectra. 
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АГ ер ате required power and cross-spectra for 
a block have been computed and stored, the sequence is repeated 
until the number of blocks requested have been analyzed or 
until an end-of-file mark is sensed on the coefficient tape. 
Then the program averages the spectral estimates both 
ET the number of blocks processed and over the analysis 
bandwidth requested in the input dala cards. (пе ваша ко 
deviation from each mean 15 computed in a similar procedure 
om the formula 
NE p 
N-1 


Б 152 
| ШОН, - R) 
а = 
EET s: the number of samples used. Other useful information 
ENSonputed and printed in either tabular form or graphic form 


as Specified bY ChE Investigator. 


Р. IBM 360/67 TAPE OPERATIONS 
1. Job Control Language 
Job Control Language (JCL) cards provide the computer 
wach information on tape mount number, tape identification, 
As position of tape; identity of tape data file to process and 
rer Of data on the tape. A typical JCL card Used соси 
meocessing would бе: 

//GO.FTO4FOOL DD UNIT=2400,VOL=SER=NPS185,LABEL=(1,SL), 
DSNAME=MCKEO1,DISP=(NEW,KEEP), DCB = (DEN=2 | 
RECFM=VS, BLKSIZE = 8204) 

GO - This group indicates data is going to be classified 


under the subsequent identify parameters. 


66 





FTO4 - This group indicates the tape is to De phoc omon 
logical unit (in this case a tape mount) four. Also the two 
digit number (04) must match the number specified in the READ 
or WRITE instruction in the program (i.e. READ (4) KMAX, NCHAN, 
DATA). The JCL does not determine whether "reading" or 
Miting" is to take place. The program statement which in- 
dicates the logical unit (tape mount) also indicates the pro- 
cess (READ or WRITE) desired. 

Eso — This group indicates the sequential number of 
passes through the tape which has occurred up to this Pointi 
In this case it is the first pass. 

UNIT = 2400 - This group indicates the particular unit 
Will be a nine-track tape. The designation for a seven-track 
tape is UNIT = 2400 -1. 

VOL = SER = NPS 185 - This group indicates that the tape 
has the external computer center marking NPS185. If a seven- 
track tape were used it could have a user specified name 
(i.e. MCKEO1). 

LABEL = MCKEO1 - This group specified the Data Set Name 
(DSNAME) of the data specified by the LABEL group. This name 
May be changed for each different file of data; however, such 
 есесйпге іс quite time consuming. It is much quicker to 
use the same DSNAME for.a whole tape. It is important that 
Een ever a data set is created under à particular name, it 
must always be referred to by the same name when reading 


from the tape. 
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DISP = (NEW, KEEP) - This group specified the disposition 
of the tape when processing (either "reading" oe "sor: dob ieee 
This implies it is a newly created tape (NEW) and it is to be 
saved (KEEP). If an old tape were to be read and the data 
saved, the selection would be DISP = (OLD, KEEP). 

DCB E = 2, ВЕСЕМ = VS, BLKSIZE = 8204; This group is 
the Data Control Block. The most often used density (DEN) is 
556 bytes per inch which is denoted by 2. This corresponds 
ШЕП the manual setting of the DENSITY dial on) the wi aces & 

ШИЕ seven-track tape drive unit. The tape record format 
(RECFM) is variable spaced (VS) when recording and digitizing 
on the seven-track tape. The block size (BLKSIZE) specifies 
the exact number of bytes in one block of data. If 2048 words 
mmo lock were considered, four bytes per word would result 

in 8192 bytes per block. If, as in the CONVERT program, two 
words from KMAX and one word from NCHAN give an additional 12 
bytes, the total is 8204 bytes per block. 

2. Multi-file Tape Operations 

Most frequently, a digital tape containing several 
files of data was generated. Under normal IBM 360 се пе 
(Пп the CONVERT PROGRAM, the JCL cards are used to indicate 
which files of data are to be processed. But, in the FIOR and 
SCOR programs, JCL cards had to be included for every file, | 
En though the file was not processed. This variation to 
the usual JCL procedures was neces. because of шаси = 


MANTE subroutine in both programs which caused the files that 


meee NOt specified for analysis to be skipped: 
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An example of the JCL cards for the normal proceso. 
of three files of data; where it was desired to skip the 
meeond file ona tape, was given by Jones [Ref. 5, pg. 46]. 

3. Multi-Volume Tape Operations 

If one long file of data extends onto another tae 
Bad processing of this long data file is desired, multizvolkme 
tape procedure becomes quite complicated when using multi- 
tape Mm its use should be avoided. The average seven- 
and nine-track tapes used can hold more than 1400 records of 
2048 samples per record. This gives a capacity for more than 


two and a half million samples per tape. 


ЕН РАЕРАКАТТОМ OF CARDS AND TAPES FOR PSD ANALYSIS 

ЖЕР UCL Cards For FTOR 

Ишет еспе ог JCL techniques for using еко 

program are taken from Jones (1971). The following Es 
assumes also that NPS 216 will be used with its KMAX set to 
256. If changes are made to KMAX and the program is taken 
from another tape (as discussed under Modifications to FTOR), 
the same argument would be valid. 


The FTOR program is stored on NPS 216 in File 1. 


The program is also storred on disk (see Wilson, et.al., 


(1969)). To use FTOR directly from the tape, the following 
Cards should be used. 
===-JOB CARD ee ewe eS کے کے کو سے سے‎ ЕЕ TO AA NT E 


//FORT. SYSIN DD UNIT = 2400, VOL = SER = NPS 216, DISP = OLD, 


// LABEL = (1,SL),-DSN = UBCFTOR 
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---JCL CARDS FOR INPUT & OUTPUT TAPES ----------------------- 

//GO. SYSIN DD* 5 

E= CONTROL CARDS FOR INPUT & OUTPUT PARAMETERS ------------.- - 
(BLANK CARD) 

E See Wilson, et. al. [Ref. T, Pe 323a mar ee 


If the FTOR deck is used, the following cards should 


EN ORT. SYSIN DD* 
----- FTOR PROGRAM DECK -------------------------------------- 
----- JCL CARDS FOR INPUT & OUTPUT TAPES --------------------- 
EGO. SYSIN DD* 
---CONTROL CARDS FOR INPUT PARAMETERS ----------------------- 
| (BLANK CARD) 
/* | 
Jones [Ref. 5 p. 70] gives à good example of how the 

input cards for a complete FTOR run should appear. His format 
 Ешпезс that the FTOR program is being input from a deck of 
Pards. 

EN c сак bration Card referred о Ваз а ве TOL 
primary and secondary channel. When the complete signal comes 
from only one "primary" input, the primary channel calibration 
Es the only one used. If an instrument, such as the 
Й іс anemometer is used, signals need to be added vectorially 
to give the actual signal. The FTOR program allows for the 
input of two separate digital signals. The calibration cards 
Permits these signals to be added vectorially and the 


eeectrum computed for a single signal. The alphameric 
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section of this card allows for signal, parameter acento 
weh as sampling rate, filter setting, analos taperehamnel 
member, etc. to be storred with the Fourier cocificicnt o mo 
 рапегіс section is printed as the title on all output 
graphs. 

onen mieht arisen ron Lhe use of this alphameric 
section when this title is applied to long data lengths. When 
OPT sections of the signal are analyzed with SCOR, this 
ШЫ е information is printed Om all сари Thus r ЕЕЕ 
bookkeeping problem occurs when thirty graphs are printed with 
EN ne title. Fortunately, the computer center delivers 
wn in a continuous roll. Identification is made by find- 
ing the starting point and then numbering them sequentially 
from ШІ. БОЛО, 

As mentioned previously, the SKIPFILE subroutine 1n 
FTOR made it necessary to include JCL cards for every file 
Exche nine-track tape. This is a variation to the normal 
IBM 360 procedure in which the JCL card specifies which file 
process. The SKIPFILE routine in FTOR, SCOR Апа PE 
fellows for the internal selection of data files, based on the 
parameters on the data cards which follow the tape JCL cards 


I the input deck. 


eee Modifications of PICK. Program 
Jones (1971) found it was necessary to make several 
changes so that FTOR would be compatible with the block size 
specified by KMAX. It was discovered that the FTOR program 


on NPS 216 was written for a block size of 256 words. The 
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change involved setting KMAX, Card 18, equal to the number of 
words in one block of data. Also the KMAX dimension statement, 
сага 23, was similarly changed. 

This change was made by inserting the new са aia 
the program deck. Then, in order that the whole FTOR deck 
need not be input each time, the modified program was written 
into File 1 of tape NPS 223. In order to reduce the number 
different tapes used, the SCOR program was written into 
File 2 of NPS 223. If too many different tapes are used, 
confusion arises insmountinge tapes. Thus, the modified FTOR 
and SCOR could be called from the same Ере оу Ето 
several cards rather than manipulating an entire deck of cards. 

fee JCL For SCOR 
. The ecards required for the "running" of the SCOR pro- 
em follow the same order as those required for FTOR. Тһе 
Nor changes are: 

LABEL = (2,SL) The program is stored in File 2 

DSN = UBCSCOR The data set name is UBCSCOR 

JCL Cards for tapes must be changed to read the 

FTOR coefficients from the FTOR generated tape. 

Control cards for analysis must be appropriate for 

SCOR program. (see Wilson, et.al. (1969) for 

details). | 

Jones, LRef. 5 p. Tll, gives a MN EC of how 
E input cards for a complete SCOR run should appear. His 
format assumes the SCOR program is input from NPS216. This 
tape can be used without the KMAX changes needed for FTOR, 
because the input to the SCOR program was the Fourier co- 


efficient tape which had a uniform format, regardless of the 


block size used in the nine-track digital input tape. 
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The selection of a uniform set or make S 
anticipated so спас сарое outputs can be overlaid orean 
comparison of spectra. If logarithmic bandwidths (exponential) 
are selected, a constant bandwidth bar appears on the log-log 
Mt due to the exponential nature oir ае асс 

Jones [Ref. 51. points out that the SUBSEQUENT ICMAX 
CARDS referred to by Wilson, et. al. [Ref. 7, pg. 36] in- 
Encre what type of graphic output is desired: spectra, cross= 
шеста ог both, The spectrum of a single channel is specified 
as if it is a cross-specturm of that channel with itself. For 
two channels (1 and 2), a spectrum of channel one and a 
Ecccprum of channel two are obtained by the following entries 
these data cards: 


First card: 1 | 2 
second card: 2 2 


ен срог5-5рестрит 15 computed only is each channel appears 
ЕГЕС list on the other channels card as seen with the first 
Para above. 
ЖОШО, Гог ЕОРГТ 

its. program is used to get a Cadvicomp plot of une 
amplitude of the Fourier coefficients as a‘ function of the 
frequency. The cards required for "running" the program 
follow the same order as required by FTOR and SCOR. The major 
pim ges required are; | 

a. LABEL = GDL) Tiis pio о оне с in File Be 

b. DSN = UBCFCPLT The data set name is UBCFCELTS 

c. JCL cards for tapes may be the same as those used 


for SCOR if the same data files are to be used. 
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а. Control cards for analysis parameters required 
mst be appropriate for FCPLT as specified in Wilson} cua 
Mier. 7, pg- 59). 

Jones [Ref. 5, pg. 72] gives an excellent example of 
how the input cards for a complete SCOR run should appear. 
Figure 26 was prepared to further clarify the JCL cards needed 
for each program. 

о. Suggestions ГОР Efficient Tape Processing 

Since the techniques involved in running this program 
бре Quite time-consuming, much attention was РТО to 
the problem of optimizing the number of runs per day. Another 
problem was the low priority for these programs. (Priorities 
Pence ЕСІП the highest priority, class Ato the lovesti telas: 
K). The low priority resulted because of two factors. кгз 
the program computation time (specified on the green JOB card 
by TIME = MM,SS, where M is minutes and S is seconds), in most 
eases, took more than four minutes of Control Processor Unit 
(CPU) time. The second factor was that the program required 
using tape input and output units which are automatically put 
inco a lower class because of the time required to mount ІШЕДІ 
and access the information on them. If the programs and the 
ése quantity of digitized signals could be completely 
stored on magnetic disks, higher priorivies could ре achieved. 
This procedure would save time if repeated analyses were to 
be conducted on a data set which did not change. In this 
study many different signals were studied, which implied™ tne 


data set varied from one signal to another. 
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c 
DU ты. DES 


//G0.FTO2F001 DD UNI T=2400-1 ,VOL=SER="CKEN4,LOnEi= 44,00 ° 
ruc, // 015Р-(040,КЕЕР),0СВ-(ПЕл-4,ӨЕСҒР-Ғ, 515126524222 
а, DD UNIT=2400,VOL=SER=NPS1P5,LAREL=(9,SL0, 





CONVERT 


===. === Additional JCL cards included for each 
file for which conversion is desirrd, 


9-Track Digitized Data Tape 





жае = < CONMERT  PROGRAMIODECK SE OSS 


DSNAME=MCKEO1,DISP=(NEW,KEEP),DCB= (DIN= 2, RELFWE US ELWS12i- 


------------- 1 OR ЕВ НОВА ЕЕХ- yT 
//GO.FT)BF000 DD UNTT2240n,UDLEZSERs «C518 5,01 35PSOLD, LABELS (1,5L), 


DSNz FCKE01,DCB2 (DENz2,RECF'- /5, BLYSTZE- 204) 


3 см | 7-------- Additional JCL cerds for sa2" File on МРОТ 2 
Fut Because of SKIPFILE routine, ‘ZL card for each 
file on wPS185 must be included, 
гтра //GO.FTO3F001 DD UNT T=2400,VOL=SE9=4F3227, 07 P= (ACW, KEEP) LABEL=(41,5L), 
77 OSN=COEF,DC8=(DEN=2,RECF "243, "3lZe= ._4) 
--------- Additional JCL for each file to De өгїї}ел-----------—------- 
on №Р5223, 


//GO.SYSIN DD + 
--------- Four CONTROL CARDS for each file for urjicr 
Fourier coefficients are to be computed, 
BLANK CARD TERMINATES COMPUTATIONS ) 


9-Track Fourier Coefficient Tape 










e DSN=UBCSCOR 

FILE 

DSN=COEF,DCB8=(DEN=2,RECFM=YS,ELKSIZE=8 204) 
on NPS185,Because of SKIPFILE routine 
JCL card for every file on NPS185S 


, muet be included. 
//G0.SYSIN OD * - 


which spectrum is desired. 
( BLANK CARD TERMINATES COMPUTATIONS ) 


insi 


PSD Graphs PSD Values 


Figure 26. Program Sequencing and JCL Cards Needed for PSD 
Analysis 


T5 


//FORT.SYSIN DD UNIT=2400,VOL=SER=NPS216,DISP=0LD,LADEL=(2,S5L), 
//G0.FTO3F001 DD UNIT-2400,VO0LESERsNPS185,DISP- (OLD, KEEP), LABEL (1,2: 


----------- Additional ӘСЕ Cerass for each Mle о 


ee POWER SPECTRUM PLOT OF ТАРЕ 001 J.D.MCKENDORIZK..OCEANO.. 
----------- Three Cards minimum for each Ғі1е Ғог-------------.-.- 





a. Program Submission Submission 

Ihe best time for processing digre ISE P оте 
debugging particular problems was found to be during the summer 
and winter break when the work load at the computing center 
ШІБ тегу light. During that time, generally, about ven runs 
mem day, depending on the length of the JOB, could be achieved: 
If all was going well only one run per day was often adequate; 
however, when problems developed, an immediate debugging run 
was essential. It was found that oversights on the part of 
 Ггосгаптег, іп the form of errors in JCL and control cards, 
the "dumb" computer was only doing what it had been instructed 
Meme, However, occassional computer malfunctions did occur. 

The next best time for program processing was on 
weekends at the beginning of the quarter. As the quarter pro- 
pressed, the weekend would offer several Е day. Toward 
вена ог the quarter, due to heavy work load, only one run: 
 агоцппа 3 A.M. resulted, if the job had been input prior to 
10 A.M. of the previous day. In general, turn around averaged 
24 hours. 

Шие most CPU used for any GrOECTAM Е  u 
ten minutes. Due to the fact that each analysis uses different 
parameters, the only rule of thumb which was found useful was 
that the FTOR procedure took about 22.5 seconds per record (2048 
samples per record) to compute the Fourier coefficients. Thus, 
400 records required 490 seconds of CPU time. The FTOR pro- 
gram was the most time consuming, thus this figure (22.5 seconds) 
can be used as an upper limit in estimating CPU time for 


Блек programs. 
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b. Stacking Programs 

Another technique used to affect faster PSD 
analysis involved submitting CONVERT, FTOR and SCOR under one 
JOB card. The reason the programs were "stacked" was to per- 
mit CONVERT, FTOR and SCOR programs to run sequentially. 
Btherwise, it was necessary to get the results of each program 
ШО шге the next program could be марис. шито mistakes were 
made in the JCL cards and if no problems existed with the 
tapes, the complete analysis would be madé In onc ттин 
problem in any step was encountered, Corrective actionia: 
meen and the remaining programs were run individually. The 
Seton as to when to stack the programs and when to run them 
individually varied with the number of files on the tapes. 
iiemiegher the number of files, the more JCL and control 


required, and the greater the chance for errors. 
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IV. " EXPERIMENTAL PROCEDURE 


The initial plan for the investigation of the noise pro- 
blen, involved digitization and spectral-analysis eng "р 
Signals. These were Eo Be Sine waves, square Waves rampo, 
ema random noise. Since these input-signal characteristics 
were readily known. the Fourier-transform and, the power 
sectra were known functions. The first signal, a l0HZ sine 
Wave, did not give the expected Spectral values. The digital 
program was checked to see if its calculations were тия 
then the A/D procedure was checked. For purposes of clarity; 
Ae naronological discussion of the experimental procedure will 


follow. 


A. ANALYSIS OF PURE SIGNALS 

ouch the first series of 10HZ sine waves showed an 
pected energy peak at 10HZ, the exponential decrease of 
energy with increasing frequency was not expected for a sine 
function. An assessment was made that the problem could be 
either in the actual A/D step or in the spectral-analysis 
programs. Before definite conclusions concerning this pro- 
blem and the noise sources could be made, it was necessary 
to gather baseline information on computer analysis of 
Mmacoretically pure signals. 

ES Computer Generated Digital Sine Function 

The program listed in Appendix Т was used to generate 


a simulated digital sine signal and the computed samples 
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were stored on a digital nine-track tape. The format of 
the tape was made to be compatable for input of this data to 
the FTOR program. Іп the sine generation program, the peak 
voltage could be yaried by changing the sine amplitude, and 
the the sampling rate could be changed by varying At. A 
standard block size of 2048 samples per block was maintained 
Bhroushout this study. The CONVERT procedure was not needed 
because the data was already in a hexidecimal format. 
E. Computer-Generated Digital, Random Signal 

Шене не SLED Was TO test the FTOR and SCOR speccral 
density analysis of a signal with a wide range of frequencies. 
To do this, a Gaussian signal, wnich has a flat Spectral 
density function, was used. The program listed in Appendix 
II generated a simulated digital random signal by "calling" 
the computer sub-routing RANDU. The peak voltage values 
Кре changed by altering the constant multiplier of YFL; 
the sampling rate could be changed by altering the sampling 
mace specified on the FTOR input data deck. The peak amplitude 
ST maintained at 1O volts; however, the different sampling 
rates were investigated. Since the data was stored on a nine- 
track tape with a format compatible with FTOR, the CONVERT 


program was not used. 


B. A/D CONVERSION AND PSD OF LABORATORY SIGNALS 

Once characteristics had been established for computer 
processing of pure control signals, the next step was to 
ENCTUlze actual signals. It was decided hai random or 


Gaussian signal would give all the information required, and 
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chus, the sine, ramp, and square waye signals сопіс Бар E 
passed. The random signal was expected to give optimum power- 


Apeetral density information for purposes of the suudy. 


l. Randon Signal 
a. Single Channel Digitization 

An Elgenco model 603 A, Gaussian noise generator 
Was used to give a random noise output. ICS charac teri otie: 
are listed in Appendix TN А frequency setting of 5Hz to 
20 KHz, attenuation schale X1.0, output voltage reading 2.62 
Vrms was input to a Khron-Hite filter model 3321 set at 2000 
HERO LOW pass max. flat and Odb gain. The filter output was 
input to the Hybrid computer through an operational amplifier 
with a ten volt gain. A sampling rate of 5000 SPS was selected. 
A total of 41 seconds of the signal was digitized onto a 
seven-track tape. An End-of-File mark was written onto the 
tape by typing the EOF option on the teletype keyboard. 

Be Dual Channel Digitization 

me second file on Che Cape mentioned above Was 
altea with data from two input channels which were digitized 
КИШ taneously. The Elgenco noise generator described above 
was used as one input, and the other input was from the random- 
noise generator which is built into the Ci 5000. 

The noise-generator was utilized with the same 
settings as above; however, the filter was reset to 1880Hz. 
Тһе random-noise generator from the Ci 5000 was Input to а 
Krohn-Hite filter also set at 1860 Hz. The outpuy ош то о 


filter was fed into a different operational amplifier with 
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a gain of 10. The sampling rate was lowered to 4000 Sps. 
ENGotal of 25.5 seconds of signal was digitized. Though ithe 
mele and dual channel cases both considered 100 records: 
Mots shorter digitizing time occurred because, Two sates 
were béing taken, simultaneously, at the rate of 4000 SPS: 


2048 Samp x 100 Blks 


a Bbk = l 
000 5апр 2. 85.5 See. 
= Бин 
сес 


Ап End-of-File mark was again written on the tape to end this 
file of data. 

The CONVERT program as used to convert from seven- 
to nine-track een The program SCOR allowed the computation 
end plotting of the spectrum of each channel and the eo- 
spectrum and the quad-spectrum of one channel with another. 

2. Random Signal and Sine Signal Е 

To determine whether a sine wave could be picked out 
of the random noise, a 1000 Hz sine and random signal were 
ес лед. An attempt to digitize a single channel of the 
Sine combined with the random signal failed due to a faulty 
patch on the analog board. PSD values were obtained, seemingly 
because the open amplifier actually picked up stray signal. The 
Sine amplitude was increased and tine second file was 
digitized. The signals in these two files were digitized at 
5000 SPS and filtered at 2000 Hz. 

Dual channel digitized samples of the 1000 Hz sine 
with amplitude +20 yolts and Gaussian signals filled the 
Third file. The ашр was increased to +30 volts and 


the two separate sine and Gassian signals digitized into the 
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meurth file. The signals in these с врана ошон ы 


were digitized at 4000 SPS and filtered at 2000 HZ. 


C. DATA FROM GEOPHYSICAL SIGNALS 

Atomspheric-temerature and velocity signals have pre- 
viously been recorded on one-inch magnetic tape by Boston 
[Ref. 1]. These signals were used as a final check on the 
System to determine if correct spectral values could be 
еп е“ед. 

The signals were reproduced on a Sangano Model 3562 FM 
tape recorder at 60 ips. The tape playback output was filtered 
SE 1000 Hz for the temperature sional, vend ак UD MIS 
velocity, the differentiated velocity and temperature signals. 
The sampling rate for the temperature-signal digitization 
was 2000 SPS and the other three signals were sampled at 4000 


DPS. The filter setting was low-pass max, flat, Odb gain. 
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V. ANALYSIS OF RESULTS 


A. PSD OF COMPUTER GENERATED SIGNALS 

ЕС Sine Wave 

Figure 2/ was the spectral plot of a computer- 

generated sine wave. The PSD values were computed for 24 re- 
meds Ot Signal giving a total signal length of 11.9 seconds. 
The total integrated power was .499 ve: However, the power in 
the 8.12 Hz band centered about 7.11 Hz had a total of .495 vV? 
(8.13 Hz X 6.09 X 107^v?/Hz). Essentially, all the power was 
шта пеа in the band between 2.05 Hz and 11.17 Hz. Though 
ШЕШЕ арреагед to be a wide band, the whole region from ll.7 Hz 
to 263.1 Hz had less than .8 percent of the total power: 


499 - .495 
. 599 


5 100 NS 

Another test run with a 200 Hz sine wave with a 
peak-to-peak amplitude of 10 volts, proved inconclusive due 
vto an error in specifying the sampling rate on the FTOR data- 
card irput. A sampling rate of 400 SPS was specified, rather 
ШІСІ 200 SPS, which was the actual rate used. The expected 
total power value of 50.у° was achiemed: however the error 
with the sampling rate shifted the frequency peak from 200 Hz 
[95:559 Hz. Though the error produced erroneous results, the 
effect of not specifying the correct sampling rate was observed. 

2% Gaussian Noise 
Figure 28 was the PSD plot of the computer-generated 


random signal. The digital samples of the signal simulated 
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a random signal of 9.8 seconds in length, sampled at 5000 SPS. 
The spectral level was found to be very flat with increasing 
meequency. Average spectral density of 3.30 x 1073 Va Не 


3 y^/Hz to a low value 


varied from a high value of 3.45 X 10 
ER 3.25 X 1073 y^/Hz. Actual variance was very low and quite 
uniform. No frequency spikes were observed and the conclusion 
was drawn that the random number generating sub-routine RANDU 
produced a true random series for at least the first 50,000 
numbers. It was noted that the prandom-number generating 


539 


рарасісу of this sub-routine is numbers before the series 
repeats itself. Since only 211 numbers were used, the full 
potential of the number generator was not fully tested. 

Figure 29 is the SPD plot of the same random series 
Eumpted at 1000 SPS rather than 5000 SPS. This changed the 
cord length to 2.05 seconds, and since the results for 
U TEecords was computed, the total length of signal sampled 
was 205 seconds. Although the sampling rate-change produced a 
higher PSD value, the mean showed little variation. The mean 
level was about 1.57 X m V^/Hz with a low of 1.62 X 107^ 
V^/Hz and a high of 1.72 JER. The variance around the mean 
evel Was very small. As expected, averaging over a fewer 
number of records (24 versus 100), the variance was higher. 
B. PSD LABORATORY SIGNALS 

i Single Channel ше 
чы orgnal Leakage into Open Ani Err ter 
Figure 30 was the PSD plot of the spectrum of 


sapnals with peak voltages of 22.0 and 23.0.) Ine та 5 
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were picked up by an open input amplifier, which was next to 
the amplifier into which the Signals were actually input: БЕ 
ореп amplifier, whose output was being digitized, acted like 
an antenna in picking up these stray signals. The plots show 
considerable consistency and several conclusions can be made. 

| The spectral peak of the +3 yolt Signal was higher 
than that of the £2 volt signals. The peak PSD was 2.68 X КЖ. 
DE Hz for the lower and 6.58 Х 1072 ne for the higher signal. 
These values multiplied by the band width, 78.1 Hz in both 
cases, gave ронег of 2.09 X 1072 v^ ang 5.12 X 1072 y? 
respectively. The power ratio produced by this difference in 
output voltage was 


28 -5 
Power Ratio = e D E DN 


NOS 
since this was a power ratio, the voltage ratio is the square 
root of the power racion, or 

Voltage Ratio = 2.45 = 1.56 

The actual voltage input into the open amplifier 
was computed from the power of each signal. The lower power 
2.09 X 10732 resulted from an input voltage of +4.5 X Lo 
Emd 5.12 X 10732 resulted from an input of +2 xX Do o 


The expected voltage ratio was computed using the 


observed input voltage values: 


Voltage Ratio = 3.0 _ 15 
en | 


The power ratio was: 


Power Ratio = 15) Ec 
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Ihe ratio of the signal уо баве Ші ИШЕ 
mine to the actual computed signal voltage gaye the actual 
Bercent of signal picked up by the open amplifier. 


20208-12 
l6 SIO Х 100 





2.3% for the 2V signal 


7.2 x 107% уз, 
3.0 100 


218 Төр Che V sona 


Assuming the leakage was coming from the voltages 
ШІ спе input lead to the OBEN amplifier Che е ае а 
Eupnal leakage would be the ratio of the voltage in the input 
lead to the signal voltage computed from the PSD. For the 
t2 V and £3 V signals the percent of signal picked up by the 
open amplifier was 2.3 percent and 2.! percent respectively: 


4.6 X 107? 


50 Х 100 


2.3% 


ED 
T.2 X 10 " x 100 


3.0 2.4% 


it the signal was leaking from Спе closed) ampli lem 
the leakage was .2 percent for both cases. 


_2 
8.6 X 10^ у 100 
5.0 


‚2% 


-2 
[2 X 10 Х 100 


3.0 { 


Thus, only a small signal leakage was observed and 
it was independent of signal amplitude. 
b. Effect of Increasing Signal Amplitude 
Figure 31 із а РӘН рТоОрФБӘНӘШІПЕ ше егесті тіне 
creasing the voltage of the signal going into the data taking 
Enmpurjer. The input ME had peak-to-peak voltages of 


БОО УО апа +30 volts. 
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The PSD for the peaks 3.92 у^ апа 9.29 V/Hz 
were computed for a bandwith of 68.4 Hz. The power was 
268V? ana 635V? respectively. This gave a power ratio of: 


Power Ratio = 9.29 
Ecc Mo cT 


и “ој саге ratio is the square root of the power ratio or 

Voltage Ratio =N 2.37 = 1.58 

Using the power Spectral density co compuertas 
voltage ratio, 268V° implied an input of 16.4 Vans апа 635 ve 
реа ап input of 25.2 S This implied peak-to peak 
meltagzes of +23.2V and +35.6V. Due to the inaccuracy involved 
in reading peak-to-peak voltages from the oscilloscope on the 
03 5000, the observed inputs of +20V and +30V could have been 
ЕО in error. 

Assuming, the inputs were of +20 om TOO ENE 
expected voltage ratio would have been: 


Кое 
20 


and the expected power ratio would have been: 

Power Ration = (1.50)* o 

The observed and computed power ratios compared 
favorable, and the difference between ovserved and computed 
peak-to-peak voltages ( 20V Vs. +23.2V and +35.6¥) were with- 
mm acceptable limits. 

The theoretical power for sine waves of +20Y and 


+30V was found from the formula; 


ve 
Ра 
e 


је 





where V is peak-to-peak voltage. This gave power of 200 y? 
and 50ү“ for the two voltage signals respectively. The 
difference between expected power level and the power level 
derived from the spectral plots was assumed to be due to the 
Еос іп reading the input signal amplitudes. 
2. Single name. Ец 

Figure 32 was the spectral plot of 40.06 seconds of 
a random signal sampled at 5,000 SPS. The spectrum level 
Was very flat to about 1.5 KHZ. Beyond 1.5 KHz, rapid decrease 
ШО Пе power spectral density with increasing frequency was 
о expected. The Зар down point occurred at 2.0 KHZ. The 
Pope of the filter, as specified in the equipment characteristics, 
was -48db/octave. The observed slope was very close to 
-96db/octave. This value would be expected if two filters had 
en cascaded, but this was not the case. The spectrum was 
. quite free of noise spikes and had no 60 Hz harmonics present. 
If 60 Hz noise was present, its level was well below -17.5db/Hz. 

The spectral level of 1.73 X 10 ^y?/Hz compared 
Mavorably with the input signal. Specifications for the 
Elgenco noise generator gave a spectral density of approximately 
DUX 10 vy, Hz at IV mas The input signal had a meter read- 


ШАР! ОҒ к The gain factor was 10. The computed spectral 


= 
density was 2.69 X 10 *y?/uz. 


-3 
DX TON N? 
Hz ) (2.6)2(10)2 = 2.69 x 107 у2/ни 


Since no accurate spectral density information on the noise 
generation was available, these values are considered to' 


compare tavorably . 
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3. Тимо Channel SD OT Gaussian Signals 
Figure 33 is the PSD plot of one of two channels of 


Fandom noise which were digitized simultaneously e Sci an 


resulted from the Elgenco Gaussian noise generator and the 


other spectrum resulted from the high frequency random-noise 
generator of the Hybrid computer. The flat spectrum of tne 
Elgenco generator was contrasted by the very а spectrum 
of the Ci 5000 random-noise 2... whose low frequencies 
contain much power. Because of its deviation from the Gaussian 
characteristics, future use of the Hybrid noise-generator 


should be avoided when a Gaussian generator is wanted. 


Б. PSD OF TURBULENCE SIGNALS 
General Signal Characteristics Found; Comparison 
With Previous Results 
a. Temperature Signal 
Figure 34 showed the PSD of the temperature signal 
recorded by Boston [Ref.1 |. A comparison was made between the 
PSD valued obtained from 56 seconds of signal that was 
digitized and analysed at the Naval Postgraduate School Facility, 
and PSD the values obtained by Boston from the same section 
e Signal that was digitized and analyzed at the University of 
meicisn Columbia facility. 
The general PSD characteristics from both analysis 
compared favorably in slope magnitude and 60 Hz harmonic peaks. 
The strong -5/3 slope region was evident on both spectra be- 


Кееп 10 Hz and 120 Hz. Significant harmonic levels were 


J5 
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Spectral Plot of Ci 5000 random Noise Generator 
output. Sampling Rate = 4000 Samp/Sec. 
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evident in the Naval Postgraduate School study at the Vives 
(60Hz) and third (180Hz) harmonics and a less significant 

ninth (540 Hz) harmonic. Boston's results showed a significant 
Nr harmonic. The third harmonic, though detectable. snowed 

EE 0minop leyel; the ninth harmonic Was nev песе е н 
Slight increase in the slope to approximately -7/3 was ob- 
served between 200 Hz and 400 Hz (600 Hz in Boston's results) 
followed by a marked decrease in the slope between 600 Hz and 
OOO Hz. 

The major difference between the two spectra is 
EEconstant power level difference. The Naval Postgraduate 
School results were higher by a constant factor of #00 which 
Bed that a voltage gain difference of 20 was present in 
mae Naval Postgraduate analysis. 7 

As can be seen from Figure 35, in which each 
Naval Postgraduate School PSD value has been reduced by- -a 
factor of 400, the spectra for the two analyses compare vary 
favorably in slope and relative magnitude. A temporal PSD 
Beets described in a later section, would have been very help- 
ful in detecting noise in Boston's results; however, time did 
Мере its production. 

b. Differentiated Temperature Signal 

A comparison Of Вов арта соо ее 
differentiated temperature signal was not undertaken in this 
study. A long section of the signal was, however, analysed. 
Figure 36 shows the cmmparison for the PSD from five minutes 


of signal obtained in this study with the PSD from three 
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different sections of data obtained by Boston (Ref. 1]. The 
meneral characterisitics were observed to agree favorable, 
though the Naval Postgraduate Sehool results were greater by 
constant factor of 300.. In Figure 37 the Naval Postgraduate 
School results were reduced by a factor of 300 and were found 
to follow closely the results from E 1'(B) and E 1'(C). The 
first harmonic of 60 Hz was found in the Naval Postgraduate 
Senool and University of British Columbia data; however, the 
surong third harmonic found in the NPS results wasn't evident 
Ші the UBC results. 
e Velocity Signal 

Figure 38 shows the PSD of a section from the 
velocity signal recorded by Boston on tape 203 (1). It com- 
pares with a section analysed by Boston and is referred to as 
U(A). The original analysis was conducted for a sample which 
was about 56 seconds in length. The original sampling rate was 
200 SPS which gave a Nyquist frequency of 1.0KHz. The Naval 
Postgraduate School analysis was conducted on the same section 
of signal, using a sampling rate of 4000 SPS and a block size 
of 2048 samples per block. Since only 56 records were analyzed, 
ре сотта! length of signal was only 51.2 seconds. Though the 
Naval Postgraduate School results are based on a slightly 
shorter signal, the PSD plots show quite similar results. 

The power level of the University of British 
Columbia analysis is lower than the Naval Postgraduate School 
analysis by a factor of 340, which implies a voltage gain 
difference of 18.4 existed between the two sets of results. 


Harmonics of 60 Hz were found in both cases in Chat cher Ins 


от 








О 20 E 


LOG FREQUENCY 


10 


Bigure 3/. Comparison of UBC and NPS Analysis of Diffrentiated 
temperature. NPS Analysis reduced by Factor of 300 
Exponential Bandwidth 


102 





: NPS Analysis of 203 (1) U(A) 
н Sample Rate = 4000 SPS 
| . Length of Sample=51.2 seconds 


-1.0 . 
E | 
o “| 
о | 
о 
o | 

КО 2.0 JE e 

о 21 

2 е 1 
су • 

4 | | 

ы, о | 

~ К © : 
> 73.0 | 25 
= о - 

> UBC Analysis © A 

= 203(1) U(A) ot | 
rS Sample Rate=2000 SBS! \ 
p: Sample length=56 sec ı Er 
= = 
m 

C5 

© 

= 


-5.0 
А 
YA 
x 
O © 
ue 
20.5 
0 1.0 2.0 3.0 


ar FREQUENCY f (Hz) 
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and third harmonic peaks stood well above the -5/3 slope of 
mae signal. A marked deviation іп results oc urn e а ЕНЕ 
PSD plots for frequencies higher than 250 Hz. The University 
er British Columbia results showed a Spectral peak at about 
eZ whereas the Naval Postgraduate School results showed 
a peak at 620 Hz. Since the Nyquist frequency of this study 
was 1000 Hz higher than Boston's results for U(A), previous 
Valles had not been reported for the frequency range 1000 to 
2000 Hz. Indications from Figure 39 are that the slope has 
decreased from the -5/3 value to a value of -2.5/3. 

2. New Results Obtained 

Based on the encouraging comparison of results obtained 
meets study with the results obtained by Boston at Universi 
SI British Columbia, the values obtained for the онегин 
шоа Of five minutes appear to offer new insight into the 
statistical properties of the geophysical processes measured. 
шәш-гайт ту, the results obtained indicate that the statistical 
maemeroies, of judiciously chosen short sections of a signal, 
Ше сіуе valid indications of the nature of over-all processes 
Пише consideration. 

With the increased capability afforded by computer 
analysis of data, new insights may be achieved toward viewing 
natural phenomena. 

a. Temporal Variations in the PSD 

A different, but not by no means new way oraiko 
playing computer PSD values of geophysical processes, is 


through temporal PSD analysis. This technique allows the 


104 


zl 0 


-2.0 
NC 
| 
N 
30 
QU 
| 
© 
Ф 
2 
QU 
=. 
O 
= 
T 
ER ЕҢ 0 
O 
[т] 
A 
2 
O 
-қ 
c5 
O 
-1 
-5.0 
mo 
0 
Figure 39. 


©, 
è 
же 
M 
ЖА! 
A т 
М 
90 
Es 
20 e EM 
LOG. o FREQUENCY 


Comparison of Slopes of Velocity PSD NPS Results 
Reduced by Factor of 340 


105 





П (1) 602 UOJSog :TRUBTS Аатоотәд рәлетлиәләріта JO STSÁTBUY SIN "Oh sandrı 


E о, 
AJNANC S43 01001 


С s u 


Ост 


E NALIATE 


- 


2098 * 
с Wo 


4 





106 





investigator the opportunity to relate Гіпсеспастспв ста ы 
Mata signal fluctuations. А section of a signal which appeared 
to have a nominal signal amplitude, as in Figure 41, had a 
low PSD value (as seen in the section of PSD plot marked 10.24 
in Figure 42. Though the characteristics of the PSD are in- 
deed determined by sampling rate, the number of samples шш егеп 
Sec ized record, filter settings, а ке том пот Бет есейе 
Пепе га] results can ђе achieved within these limitations. 
inter е епт уок ауру гелкге ст 
(digitized block of data) was maintained throughout the analysis 
И ег UBC FTOR and UBC SCOR, the PSD of sequential sections 
of the signal could be computed. Temporal variations in the 
PSD of temperature and velocity signals were drawn from the 
Pesultvs obtained by sequentially analyzing a constant number 
 несогсів. Plots were also drawn from the results of 
Кели ап increasing number of records, beginning with the 
Eu record. These plots showed that the fluctuating PSD 
Er creme more stationary when more and more data is analyzed. 
(1) Temperature. 

Figure 42 shows the temporal variation of 
the temperature PSD. A total of 300 records were analyzed to 
give a total time of 307.2 seconds. To compute the PSD values, 
thirty passes were made through the Fourier coefficients with 
the PSD being computed for each set of ten records. The time 
difference between.each PSD was 10.24 seconds. 

Figure 43 showed how the temporal variations 


in the PSD are smoothed out by taking successively longer 
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Brush Recording of Undifferentiated and Differentiated Temperature 


Figure 41. 


Signal (203(1) El and El') 





Teusts Jo spuooeS he OT UOLT TOJ 
pe4nduo3 dSd *Teusts sangeasdus] Orasydsougy FO SUOTYJETAIEN Te uaodua J, 


(2) Y AININOIYI 7901 


35:63) 13 WNYLIJAS “907 








toy DNL 
os Sitti hed - - - 
мн ' 












7 


ИЛА ЖАМУ 


А] я ы | С 
fs ЛОМЕ ИУ SÓN 








/, ) РА А ex ISX А 
eo Ж | / / N i K / 7 и“ 
NES STINK \ " МУ 7 


SEN AM | 


“ен SUNSTA 


109 





STSÁTBUY dSd доу рэзАтеиу грдоогн Jo дэашим JUTSPTIADUT JO 109333 ‘EH Sandra 


(2839 3 XA9N3nO344 7903 


0% Se ос sí ot s'i оі S: 


04- 
Е = 
О Lo 
NN џ 
к EN ES a 
a m \ 8 
| М. M ! HU ` AJ 
ођ St ot .-. : ; ee \ 
Ls NAN BS \ е : 01- E 
К" we ` ‘ \ ER ENS \_ zs 
5229 IN N : \ N N D N г) 
VAL N у ES AS N = X 
N \ р VEN SN ~ x (7 
an M zn 
\ ео ~ N E aS = * Mo S ыл 
E ү SS ~ И г 
5,4 Nu N EN B N З ка N N, ng 
| A а os E = ES AS к Sy 
Ме AM Te ~ No ES E 
NS К \ B. N \ N | | 
Кы а 2 N by i sS N i 
N ч. ч М. m Bow у oÇ 
~ " а e 


AW NN . 0v | 


110 





records of data. The first PSD was computed Proms е ои и 
starting with record number one. The second PSD along the 
time axis was computed from 20 records, starting with the 
first and so forth. The last PSD gave the average of a total 
of 102.4 seconds of temperature fluctuations, which was com- 
pured from 100 records. The smoothing абз was quite evident 
from the fact that the first PSD was a much lower value, 
enough, its effect was quickly smoothed by the averaging pro- 
cess. 

(2) velocity. 

Figure 44 shows the temporal variation of 
the loci ty PSD. A total of 300 records were analyzed. Bes 
cause of the higher sampling rate of 4000 SPS, the total 
length of signal shown ЕЕ LSI Seco. 


2048 Samp/Record 0 “ 
КЕШ лаш ы = х 300 Records 


Штерту separate PSD values were computed by sequentially 
Пас, 10 records at a time. The time difference between 
MES D curve was 5.12 seconds. 

Figure 45 showed the smoothing effect 
achieved by analysing an increasing number of records of' the 
velocity signal. The overall trend in smoothing was not as 
apparent for this signal due to the fact that the original 
section of signal analyzed gave a good statistical description 
ENSE process. 

b. PSD Five Minute morto rena 
Figure 46 showed the temperature PSD values for a 


long length of data (5 minutes) compared with data from a 
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Comparison of 56 Records and 5 Minutes of Temperature 
Signal 
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STE Section of signal. The two lengths compare very closely, 
implying that for this situation, the shorter ener вы 
minute would have given statistics representative of longer 
sections, (for the frequency range 1 Hz to 1KHz). 

Figure 47 showed the slope characteristics for 
the slong record length (5 minutes). A definite increase 


(-7/3) in the slope was noted from about 70 Hz to about 300 Hz. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Сто Томапе conclusions were reached concerning Chio 
problem of the possible input of noise into the Analog-to- 
Digital conversion procedure and into the PSD analysis pro- 
Esdure with signal inputs Of real data. 

URESD Programs 

No Significant noise sources were found to exist with- 
поште computational program of the Naval Postgraduate School 
FFT package. The PSD computations for noise -free sine and 
random signals agreed very well with theory. 

The programs FTOR, SCOR and FCPLT were found to be 
EE enu Lor obtaining power spectra from oe quantities 
of time series data. 

In the CONVERT program used in a previous study, it 
КЕШ оппа that a factor was missing which changed the octal 
base number to the hexidecimal. A corrected version of the 
program was used and PSD results indicate the procedure is 
mumecioning correctly. 

PSD values from four atmospheric turbulence signals 
which were compared with results obtained from other 
feompucational facilities showed very close correlation. 

2. Analog-to-Digital Conversion 
No significant noise sources were found vo ре inter- 


fering with the digitization procedure carried out on the 


TEF 





Hybrid Computer. Patch board noise which Had Бе M и 
reported as excessive by Jones, [Ref. 5] was noted; however, 
и ртегепсе in the PSD plots of sine andirandonk NE 

were not EET Even the noise picked пр би сре 2680 


2 Emo 


amplifier had a very low level, on the order of 1079 V 

Early in this study it was discovered that the 
digitization procedure was missing several data samples at the 
mma of each record. The problem was due to internal delays 
within the XDS 9300 which caused several data samples to be, 
omitted. The computer laboratory staff revised the digiti- 
Zation program to include a cross-check between each sample 
and lapsed time. 

it is now possible to sample at a rate of at least 
DRODZSPS aná check each block for missing data. Results 
obtained from PSD analysis indicate the problem has been 
corrected. | 

O PSD Analysis Procedures 

The IBM 360/67 was found to be quite capable of pro- 
cessing large quantities of time series data. The operational 
routine in PSD analysis has been improved. Methods for 
affecting faster "turn around" have been developed. 

Several methods Have been developed for statistically 
checking the digital values on tapes and plotting the values 
лиле Calcomp plotter. 

4. PSD Analysis of Turbulence Signals 


A definite correlation was noted between the temporal 


BSD plots of both temperature and velocity and the original 
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analog signal from which the PSD results were Com е ти 
tended to further support the conclusion that the PSD pro- 


grams were working properly. 


B. RECOMMENDATIONS FOR FUTURE WORK 

Du to the fact that this series of programs haszprer 6 
be an extremely powerful tool for signal analysis, its future 
use should be vigourously persued. The programs'use in tur- 
bulence anlysis has been well established; however, its ap- 
plication in any field which employs PSD techniques should be 
followed. 

The time-varying PSD analysis of turblence signal should 
be persued. This technique showed interesting possibilities 
for future work. Digital tapes of four turbulence signals 
were developed and can be easily assessed for future analysis 
of the signals. Though a correlation between temperature and 
velocity was noted, new digital tapes will have to be generated 
with the dual channel digitization procedure in order to get 


Cross—spectral values for these signals. 
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APPENDIX C 


Equipment Specifications 


ДЕ. 


Gaussian Noise Generator 
Manufactured by Model No. Frequency Ranges Output Spectrum 
Elgenco, Inc- 603A 5Hz-20KHz SO 
(10Hz-500KHz) 


Output Level Spectral Density mv/ Hz 
ZOKHz Setting 5 mv/ Hz 
3 V @l V output and 
E 20KH "Range 
Filter 
Manufactured by Model No. Frequency Range 
Khron-Hite, Corp 3340 O01 TG G0 MU 


Begequency Accuracy Pass Band Gain Aetenuation Siope 


227 odb or 20db -48db/octave 
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PSD Values from Computer Generated 10 Hz Sine Sine Wave 
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PSD values from Computer Generated Random Signal. Sampling Rate = 5000 SPS 
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PSD Values from Computer Generated Random Signal. Sampling Rate - 1000 SPS 
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PSD Values for Signal Leakage into Open Amplifier 
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PSD Values for Signal Leakage into Open Amplifier. Increased Signal Amplitude 
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PSD Values for a Real 1000 Hz Sine Wave. Amplitude +20 Volts. 


APPENDIX I. 
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